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ABSTRACT 

We argue that one of the basic assumptions of the Bondi accretion process, that the accreting 
object has zero pressure, might not hold in many galaxies because of the pressure exerted by 
stellar winds of star orbiting the central super massive black hole (SMBH). Hence, the Bondi 
accretion cannot be used in these cases, such as in the galaxy NGC 3115. The winds of these high- 
velocity stars are shocked to temperatures above the virial temperature of the galaxy, leading 
to the formation of a hot bubble of size ~ 0.1 — 10 pc near the center. This hot bubble can 
substantially reduce the mass accretion rate by the SMBH. If the density of the hot bubble is 
lower than that of the interstellar medium (ISM), a density- inversion layer is formed. Adding to 
other problems of the Bondi process, our results render the Bondi accretion irrelevant for AGN 
feedback in cooling flow in galaxies and small groups of galaxies and during galaxy formation. 



INTRODUCTION 



It is widely accepted that feedback powered by active galactic nuclei (AGN) has a key role in galaxy 
formation and in coo ling flows in galaxies and in clusters of galaxies. In galaxy formation AGN feedback heats 
and expels gas (e.g.. lOstriker et al.ll2010l and references therein), and by that can determine the correlation 
between the central super-massive black hole (SMBH) mass and some properties of the galaxy. In cooling 
flow clust ers jets launched by the SMBH heat the gas an d maintain a small, but non zero cooling flow (see 
review by McNamara fe Nulsenl 2007 . 2012 ; Fabian 2012 ): this is termed a moderate cooling flow. 



There is a dispute on how the accretion onto the SMBH occurs, in partic ular in cooling flows. One camp 
argues for accretion to be of h ot gas via the Bondi accretion process (e.g., Allen et al. 20061 : Russell et al 



2010c iNaravan fe Fabianll201lf ). while the other si de argues that the accre tion is of dense and cold clumps 
in what is termed the cold feedback mechanism ( Pizzolato fc Soker 2005h . The cold feedback mechanism 



has been strengthened recently by observations of cold gas and by more detailed studies ( Revaz et al. 2008; 



Pope l2009r IWilman et al. 2009. 2011; Nesvadba et~al1 


|2011 


Cavagnolo et al. 


2011; 


McCourt et al. 20121 


Sharma et al. 2012; Farage et al.l 


2012; 


Kashi et al. 


2012ft. 



•2012a. b 



The Bondi accretion process, on the other hand, suffers from two problems. The first problem is that in 
cooling-flow clusters the Bon di accretion rate is too low to account for the AGN power (e.g.. lMcNamara et al 



2011 ; Cavagnolo et aljboil ). The second is that there is no time for the feedback to work ( Soker et al.l 2009T T 
This is because the inflow time scale from even a moderate distance from the center, ~ 1 kpc, in the Bondi 
accretion process is much longer than the radiative cooling time of the gas. The gas at large distances has 
no time to communicate with the SMBH before it cools. 

In this paper we poin t out yet another problematic point with the Bondi accretion process. In a recent 
paper. IWong et al.l (|201ll) resolved the region within the Bondi accretion radius of the SO galaxy NGC 3115. 
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If the density and temperature profile is interpreted as resulting from a Bondi accretion flow onto the 
M B h = 2 x 10 9 M Q central SMBH, the derived accretion rate is M B = 2.2 x 1O _2 M yr _1 . They note that 
for a radiation power of 0.1Mb c 2 , the expected accretion luminosity is six orders of magnitude above the 
observed upper limit. They attribute this to a process where most of the inflowing gas is blown away. 

We take a different view. We argue that the Bondi accretion flow is not relevant for the conditions in 
typical galaxies and clusters of galaxies. The reason is that one cannot assume a zero pressure at the center, 
either because of stellar winds or because of jets blown by the AGN. 



2. THE PRESSURE OF STELLAR WINDS 

The pressure exerted by stellar winds of high velocity stars (i.e., moving much faster than the dispersion 
velocity in the galaxy) with an average mass loss rate per star of can be calculated in two limits, which 
basically lead to the same result. First we calculate the pressure by considering the total outward momentum 
flux at radius r. The relevant velocity in general is not that of the wind relative to the star, but rather the 
velocity of the star under the gravitational influence of the SMBH, 



Vio 9 M y v p c 

This holds as long as the SMBH gravity dominates that of the galaxy. In NGC 3115 that we study in more 
detail in section [31 for example, the SMBH gravity dominates that of the galaxy to a distance of ~ 30 pc as 
the black hole mass is Mbh = 2x lO 9 Af . Let stellar winds from high- velocity stars dominate the pressure 
inside a sphere of radius R. The pressure exerted by the wind on a surface of radius R is approximately 
given by adding the ram pressures of winds from all stars inside the sphere of radius R, 

AnR 3 1 

P m *(R) ^ n*rjm*(R)u* - - , (2) 

3 AttR z 

where is the stellar density in the center of the galaxy, and rj is the fraction of the mass lost by stars that 
is shocked and heats up. In all our expressions the stellar mass loss rate appears as rym*. 

Some of the mass lost by stars will form dense clumps that will cool rapidly even if being shocked, or 
will not even be shocked. This is particularly true as most of the mass is being lost by asymptotic giant 
branch (AGB) stars that have dense winds. The thermal pressure of the ISM in the center will cause part 
of the winds' gas to form dense clouds. The average mass loss rate is calculated as follows. A solar-like star 
loses ~ O.5M over ~ 10 10 yr. Considering an old population of stars, the mass loss rate is lower even. More 
accurately, most of the mass loss is due to AGB stars, which live for ~ 10 7 yr, and lose mass at an average 
rate of ~ 10~ 7 Af Q yr -1 . During the final stages of the AGB the evolution is faster and the mass loss rate is 
higher. If there is a young stellar population, the total mass loss rate can be much higher. The ram pressure 
will not increase much beyond few pc because the stellar density decreases. 

An alternative point of view would be to express the pressure as (roughly) the energy density of the 
shocked stellar wind. We also assume a constant pressure and density inside this sphere. This is not the 
case, but it is adequate for our goals. We can calculate the rate of energy input and multiply by the time it 
takes the hot gas to leave the inner region 

Tosc = W*' (3) 
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where (3 < 1 takes into consideration that the hot gas at the center escapes at velocity lower than the escape 
velocity. The stellar wind pressure in this case can be written as 

P =-^— (4) 

3 V fill* 

where the energy deposition rate is 

E = J ^-n*r]m*ul^ 4nr 2 dr = 2nGMBHrjm* J n*(r)rdr. (5) 
Scaling the different quantities and assuming a constant stellar density we find 

P..M = 3 x 10-,-', ( " 2 ( *) " 2 ( r -±« ) ( *= ..„_, ) org (6) 



Vio 9 M y V 1 P c / V 5 x 105 p c ~ 3 / v 10 ~ 10M o y r 

where the stellar densit y is scaled by the average stellar density within ~ 3 pc from the center of NGC 3115 
(jKormendv et al.lll996T ). Equations (j4|) and ([6]) are more accurate than equation (f2j) when the radiative 
cooling time of the colliding stellar winds is larger than the escape time T esc , which is the case here due 
to the high-temperature low-density post-shock stellar winds. Therefore, from now on we will refer to the 
hot gas region formed by the shocked stellar winds as the hot bubble. For a constant stellar density within 
radius r, we find P e * = §/3 -1 P m *. If the stellar density drops to zero at some radius r z (a nonrealistic ideal 
case), the pressure beyond r z will drop like (r/r z )~ 2 . 

The average density of the hot shocked stellar wind is given by 

47r \ 1 R f^' 

— R 3 ) 77777,*— — — / ATm*(r)r 2 dr (7) 
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The flow structure is schematically drawn in Fig. [TJ 



THE CASE OF NGC 3115 



At the Bondi radius R B ^ 210 pc of the galaxy NGC 3115 the ISM pressure is P(R B ) = 2 x 
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erg c m 3 , the electron number density is n e (Rs) — 0.02 cm , and the temperature is T(Rb) 



10 

3.5 x 10 6 K ()Wong et al.ll201lh . The Bondi radius is given by 



Rb 



2GM bh 



220 



Mi 



BH 



2 x 10 9 Mr. 



T 



3.5 x 10 6 K 



pc, 



(8) 



where c s is the sound speed in the undisturbed gas. The temperat ure and electron d ensity increase inward, 
reaching values of T 2 o — 10 7 K and n e 2o — 0.3 cm -3 at r — 20 pc (I Wong et al.ll201ll ; no values are given at 
smaller radii). We also note that in NGC 3115 the BH gravity dominates that of the galaxy to a distance of 
~ 30 pc as the black hole mass is Mbh = 2 x 10 9 M q . 



The average density and pressure of the hot bubble according to equations ([7]) and 
Fig. [2] for a SMBH mass of Mbh = 2 x 10 9 M q , and a stellar density given by 



are drawn in 



5 x 10 5 pc 



_ 3 I 1, r < 3 pc 

1 (r/3 pc)" 3 , r>3pc, 



(9) 
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Fig. 1. — A schematic drawing (not to scale) of the flow structure where a hot bubble, formed by stellar 
winds of high-velocity stars orbiting the central SMBH, exerts pressure on the ISM residing outside radius 
R. If the density in the hot bubble is lower than the ISM density, the flow at R is RT-unstable and a 
density-inversion layer is formed. 



and for (3 = 1 (eq. [3]) and 77 — 0.1 (eq. [5]). The density within r — 3 pc is from Kormendv et al. ( 19961 ). while 
at r > 3 pc is our assumption. The value of the mass loss efficiency, which is the fraction of the mass lost 
by stars that ends up as hot gas in the hot bubble, is chosen as r\ = 0.1 to more or less match the pressure 
and density of the ISM at r — 20 pc. It is a parameter of the model that should be typically in the range of 
r~> 0.1 — 1. The temperature that is calculated from the pressure is also drawn on Fig. [2] Beyond ~ 30 pc 
the average temperature is only ~ 2 as large as the virial temperature of the cluster, and our assumptions 
of a hot bubble become inadequate. 

The following conclusions emerge from Fig. [5] (1) The pressure of the shocked stellar winds of the high- 
velocity circum-SMBH stars is larger than the ISM pressure near the center, even for a mass loss efficiency of 
only T) ~ 0.1. This accounts, we argue, for the accretion rate of NGC 3115 being much lower than the Bondi 
accretion rate (jWong et al.ll2011l) . (2) At the center, r < 3 pc, the rate of mass loss into the hot gas per unit 
volume is x = (?i*?7rh*) c = 5 x 1O -6 M0 pc -3 yr . Even if this value is ten times lower, a hot bubble with 
pressure larger than the ISM pressure of NGC 3115 can still be formed. (3) For \ < 1O _5 M0 pc -3 yr -1 the 
hot bubble's density is lower than that of the ISM. This structure is Rayleigh- Taylor (RT) unstable. This 
structure is analyzed below. 
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Fig. 2. — The average density and pressure of the hot bubble according to equations ((7} and ([5]), as well 
as the temperature that is calculated from the pressure for the stellar density profile given in equation [5] 
The escape velocity parameter is j3 = 1 (eq. [3]), and the mass loss parameter of 77 = 0.1 (eq. 0) is taken to 
crudely fit the ISM properties of NGC 3115 at r = 20 pc. 



4. A TENUOUS HOT BUBBLE FORMED BY STELLAR OR AGN WINDS 



We found above that in some cases the hot bubble that formed by the stellar winds of circum-SMBH 
high-velocity stars can have a lower density than the ISM while its pressure is about equal to the ISM 
pressure Pism- This situation is prone to RT instability. The same might hold for AGN winds. The power 
of the winds that is required to form a hot bubble that can support the ISM is 

"'-..1 - -f'^Vr-J =>> 10 37 /?f ln7 ^" e .„_^ (l^-] 1/2 (Az) 3/2 ergs- 1 , (10) 



3 

2-Pism 
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where the escape time r csc is given by equation Q, and V is the volume of the hot bubble. This implies that 
even a very weak AGN wind can form such a bubble. With an efficiency of 1%, namely, W W md = O.OIA/bhc 2 , 
the required accretion rate is Mbh = 5 x 10- 8 M q yr" 1 . For comparison, we note that the Chandra upper 
limit on the luminosity of NGC 3115 is - 10 38 erg s" 1 (jPiehl fc Statle jl2008h . 



The flow structure considered in this section has the following properties. The hot bubble is continuously 
supplied by hot gas from the shocked stellar winds or the AGN wind or jets. A pressure equilibrium is 
maintained between the hot bubble and the ISM, and a structure of a hot tenuous gas supporting a denser 
and cooler gas is achieved. This structure is RT unstable. Such a structure, we claim, is similar to the de nsity 



inversion found in the outer atmosphere of red giant stars (e.g., Harpaz 1984 : Frevtag fc Hofner 20081) . but 



not identical. At the outer edge of the recombination zone of hydrogen in red giant stars the convection heat 
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transfer becomes less efficient. The requirement to transfer energy leads to a steep tempera ture gradient 
that in turn causes a density inversion, i.e., the density increases outward (e.g., Harpazl 1984 ). This occurs 
in the convective region, which is already unstable. In the density-inversion layer in stars, therefore, cold 
convective cells fall and hot convective cells buoy outward. We suggest that the same process occurs in the 
flow structure discussed here. 

There are some basic differences in the properties of the density-inversion layers of stars and of the case 
studied here. The main differences are that the hot gas in our case buoys to large distances, and fresh gas 
from stellar wind or the AGN replaces it. Also, the entire region is optically thin, unlike stars where it is 
optically thick. In stars the width of the density-inversion region is determined by heat transfer requirements, 
whereas in our case it is determined by dynamics, mixing, and local heat conduction. In stars the density 
scale height is not much shorter than the pressure scale height l p . The size of the convective cells is taken 
to be of the order of the pressure scale height. In our case the density can change by an order of magnitude 
from the inner tenuous region to the denser outer ISM, and we expect the RT instability to break the cells 
to smaller cells. We therefore take the size of the rising and falling gas elements to be R c ^ l p ~ r. 

We take the density-inversion zone to be of the order of the pressure scale height (in stars it can be 
much smaller) . For a central gravity source the pressure scale height for a constant temperature is given by 



l p = r 



a 



u*(r) 



(11) 



where Ci is he isothermal sound speed, and u* is given in equation ([TJ. For an average temperature of the 
hot gas considered here and the ISM, we can take l p ~ r. Therefore, we assume first that the width of the 
density-inversion layer is Ar^ ~ r. 

Consider then a spherical parcel of gas (a blob) of radius R c and density of pt moving with a terminal 
velocity Vj, through an external medium of density p e . The buoyancy force on the blob is 



Fb = (p e - Pc) ^Rlg- 



(12) 



The drag force on the bubble is 



F d 



C D nRlp e Vt, 



(13) 



where Cd 
is 



0.75 (|Kaiserl l2003). Assuming p c <C p e and taking g — u^/r, the terminal velocity of the bubble 

V2 / D \ 1/2 



(14) 



where in the second equality we identify the terminal velocity as the velocity by which the hot gas escape 



from the hot bubble outward, with 



/3 ~ 0.6 



Rc 

O.lr 



1/2 



(15) 



Complex processes take place in the density-inversion layer. (1) Heat conduction time scale over a 
distance of Ary = R c ~ 0.1 pc and a temperature difference of AT = 10 7 K, is few x 10 yr. This is shorter 
than the fall time of a dense clump from ^ 1 pc. Therefore, the hot bubble gas heats the clump by heat 
conduction. Closer to the center, the clump will be shredded to smaller cells. Hence, before the dense 
ISM clumps can reach the center it will be eva porated. This is not tru e for denser and cooler blobs that 
fall inward, as in the cold feedback mechanism ( Pizzolato fc Soker 20051) . (2) Because of the stellar motion 
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and/or AGN activity, the density-inversion layer is expected to be more chaotic than just a RT-unstable 
region. There will be vortices that will increase mixing, namely, reduces the effective value of Att- 



5. DISCUSSION AND SUMMARY 



We studied the pressure exerted by the winds of circum-SMBH high-velocity stars on the surrounding 
ISM. Wc found that in some cases this pressure is significant and can substantially suppress the inflow of the 



ISM r elative to what a simple Bondi accretion would give. Our result can explain the finding of IWong et al 



(|201 lh that the Bondi accretion rate calculated by them from the ISM density and temperature is six orders 



of magnitude above the observed upper limit on the accretion rate in the SO galaxy NGC 3115. 

In section|3]we quantitatively examined the situation in the galaxy NGC 3115. Shocked winds of circum- 
SMBH high-velocity stars form a bubble of hot gas whose pressure is significant, as evident from Fig. [2] 
There are some uncertainties in the model, such as the exact behavior of the stellar mass loss, trajectories 
of stars around the SMBH, and the stochastic behavior of the post-shock stellar winds. Some of these will 
be studied in future numerical simulations. However, the result that the stellar winds cannot be ignored 
is robust. For some values of the parameters we found that a situation might arise where the hot bubble's 
density is lower than the ISM density. In this case, Rayleigh- Taylor (RT) instability takes place, and a 
density- inversion layer is formed (see schematic description in Fig. [I]). Although hot tenuous gas buoys 
outward and dense ISM gas moves inward, the density-inversion layer itself continues to exist. The ISM gas 
is heated near the center and accumulated to the hot bubble. 

Our result is more general in showing that in many cases the Bondi accretion process does not work 
because one of its basic assumptions, that there is no central pressure, breaks down. This is one of several 
reasons why the Bondi accretion model may not apply in some cases (see section [T]). 

Finally, we note that our model may be relevant for active galaxies where the hot bubble might be 
formed by the AGN jets or winds. For typical values of AGN jets and winds the hot bubble density will be 
low, and a density-inversion layer will be formed. We expect this process to be of high significance in the 
process of AGN feedback acting in young galaxies. Barring Bondi-like accretion, dense and cold clumps in 
the ISM can still flow inward and feed the SMBH. Namely, AGN feedback mechanisms require the feeding 
to be by cold clumps, i.e., a cold feedback mechanism. 

This research was supported by the Asher Fund for Space Research and the E. and J. Bishop Research 
Fund at the Technion, and the Israel Science foundation. 
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